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Abstract 

We carried out a study of neutrino detection at the experimental fast reactor JOYO using a 0.76 tons gadolin- 
ium loaded liquid scintillator detector. The detector was set up on the ground level at 24.3 m from the JOYO 
reactor core of 140 MW thermal power. The measured neutrino event rate from reactor on-off comparison was 
1.1 l±1.24(stat.)±0.46(syst.) events/day. Although the statistical significance of the measurement was not enough, 
the background in such a compact detector at the ground level was studied in detail and MC simulation was found to 
describe the data well. A study for improvement of the detector for future such experiments is also shown. 
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1. Introduction 

Reactor neutrinos have been playing an important 
role since its first discovery in 1956 [1] for the progress 
of elementary particle physics and to deepen our under- 
standing of the nature. Now the reactor neutrino de- 
tection techniques have become mature after a number 
of reactor neutrino experiments so far performed 0] . 
Research and development of compact reactor neutrino 
detector utilizing the up-to-date technologies have be- 
come active recently [4] with an idea of using it as a 
monitor for Plutonium breeding in reactor cores 10] 
and as a very near detector to calibrate reactor neutrino 
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13 flux for long baseline reactor neutrino oscillation exper- 

14 iments. 

15 1.1. Reactor neutrinos 

16 In operating reactors, 235 U, 238 U, 239 Pu and 241 Pu per- 

17 form fission reaction after absorbing a neutron. The As- 
ia sion products are generally neutron-rich unstable nuclei 

19 and perform jS-decays until they become stable nuclei. 

20 One v e (anti-electron neutrino) is produced in each /3- 

21 decay. The energy of the reactor neutrinos corresponds 

22 to yS-decay energy of a few MeV. Roughly 6 v' e s are pro- 

23 duced in a Assion reaction along with ~ 200 MeV of 

24 energy release, resulting in 6 x 10 20 v' e s production per 

25 second in a 3 GW,/, power reactor. 

26 1.2. Nondestructive Plutonium Measurement 

27 Main components of reactor neutrinos come from 

28 235 U and 239 Pu Assions, and contributions of 238 U and 

29 241 Pu are much smaller than those nuclei. Along with 

30 the burn-up of the core, 235 U is consumed and 239 Pu is 
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Isotope 


v(>1.8MeV) 


q (MeV) 


Contribution 








@JOYO(%) 




1.92±0.02 


201.7±0.6 


37.1 


238 u 


2.38±0.02 


205.0+0.9 


7.3 


239p u 


1.45+0.02 


210.0±0.9 


51.3 


241 Pu 


1.83+0.02 


212.4+1.0 


4.3 



Table 1: Number of v e per fission with the energy above 1.8MeV @] 
and energy release per fission for major isotopes in nuclear reac- 
tors 0J . 



31 'breeded' from U through neutron absorption and fi- 

32 dacays. Because 239 Pu can be used for nuclear explo- 

33 sion, it is an important object of strict safeguard reg- 

34 ulations. Therefore, it is important to monitor reactor 

35 operation and track the plutonium breeding. Interna- 

36 tional Atomic Energy Agency (IAEA) watches reactors 

37 in the world with surveillance cameras, reviewing oper- 

38 ation record, etc. Because it is impossible to hide the 

39 neutrinos, it could be a powerful tool to monitor the re- 

40 actor operation, in addition to the traditional monitoring 

41 methods [@] . 

42 The reactor neutrino monitoring has a potential to 

43 non-destructively measure the plutonium amount in the 

44 core. 

45 Table[T]shows the energy releases and expected num- 

46 ber of emitted v e 's above 1.8 MeV per fission, and av- 

47 erage ratio of fission in the JOYO core for major iso- 

48 topes in nuclear reactors. As shown in the TableQ] 235 U 

49 produces significantly more neutrinos than 239 Pu. Com- 

50 bining the neutrino flux and thermal power generation, 

51 there is a possibility to measure Plutonium amount in 

52 the core. This is simply depicted by the following equa- 

53 tions assuming the fuel is made up only from 235 U and 

54 2 39 PU. 

^235-^235 + ^239^239 = Pth, (1) 
^235-^235 + V239-F239 = Ny e (2) 

55 where, 235 and 239 represent 235 U and 239 Pu. F x is 
se the fission rate of the nucleus-x in the core, q x is the 

57 energy release per fission. v x is the expected number of 

58 emitted v e 's per fission, N Ve is the total emission rate of 

59 v e . A small contribution from 238 U and 241 Pu is ignored 

60 to simplify the calculation. 

61 The fission rate of 239 Pu is calculated from those re- 

62 lations and the values of the parameters, and the 239 Pu 

63 amount in the core can be calculated from the fission 

64 rate. 



65 1.3. Compact neutrino detectors 

ee As R&D of compact neutrino detectors, an exper- 

67 imental program led by Lawrence Livemore National 

es Laboratory (LLNL) and Sandia National Laboratories 

69 (SNL) measured neutrino energy spectrum at a short 

70 distance from a 235 U-rich reactor with a thermal power 

71 of 3.4GW,/,, San Onofre Nuclear Generation Station 

72 (SONGS), and indicated feasibility of the neutrino mon- 

73 itoring |9]. On the other hand, further R&D studies of 

74 detector design and materials are still necessary to real- 

75 ize a compact detector operation above ground for prac- 

76 tical use as a reactor monitor with the neutrino detec- 

77 tion. Considering the neutrino interaction cross-section 

78 on proton target (inverse /3-decay, 0(1CT 43 ) cm 2 , see 

79 Section[3]) and compact detector size, the detector must 
so be set at a short distance (less than a few tens of meters) 
si from the reactor core to accumulate enough statistics 

82 for monitoring. In addition, feasibility of the measure- 

83 ment at ground level is required for the monitor consid- 

84 ering limited access to the reactor site, while the previ- 

85 ous measurements of neutrinos were operated at under- 
se ground to reduce cosmic-ray muon background. There- 
87 fore, the detector must be designed to be able to reduce 
sa external backgrounds, e.g. cosmic -ray muons and fast 

89 neutrons. 

90 We constructed a 0.76 tons gadolinium loaded liq- 

91 uid scintillator detector as a prototype of KASKA de- 

92 tector [10] and we reused it to take part in such R&D 

93 efforts [11]. The detector was set up at 24.3m from 

94 Joyo experimental reactor core whose thermal energy 

95 was 140 MW 11211 . Unique points of this experiment 

96 are, (1) the reactor power is much smaller compared 

97 with the ones so far used to measure the neutrinos, (2) 

98 the detector is located above ground, (3) the reactor was 

99 a fast reactor, so that the neutrinos came mainly from 

100 Plutonium. The main goal of this experiment was to 

101 distinguish reactor-on and off by neutrinos under this 

102 unfavorable conditions. One of the possible safeguard 

103 applications is to monitor small reactors to prevent them 

104 to be hiddenly operated to make plutonium. The points 

105 (1) and (2) of this experiment are useful to study such 

106 a possibility. As for (3), neutrinos from 235 U-rich light 

107 water reactors have been measured lfl3lllfl4ll . while ob- 

108 servation of neutrinos from 239 Pu-rich fast reactor has 

109 not been reported yet and this experiment could have 

110 been the first detection of the fast reactor neutrinos. If 

1 1 1 energy spectrum of fast reactor neutrinos is measured in 

112 the future, V235 and V239 can be determined separately by 

113 comparing the 239 Pu-rich neutrinos and 235 U-rich neu- 

114 trinos. This experiment is a good practice to perform an 

115 experiment at a larger fast reactor in the future to mea- 

116 sure the Plutonium-rich neutrino spectrum. 



2 



1600 
1400 
1200 
1000 



06/10/03 06/11/02 06/12/02 



Figure 1 : Fission rate of each fissile element as a function of time from 
the 4th to 6th operational cycles of experimental fast reactor JOYO. 
Four lines correspond to 239 Pu, 235 U, 238 U and 241 Pu as indicated in 
the figure. A period used for data analysis is also shown. 



2. Experimental fast reactor JOYO 

The experimental fast reactor JOYO, whose thermal 
power is 140 MW, is located in Japan Atomic Energy 
Agency (JAEA) Oarai Research and Development Cen- 
ter in Ibaraki prefecture, Japan. The JOYO reactor is a 
sodium-cooled fast reactor built as an experimental re- 
actor to promote commercialization of fast breeder re- 
actor development 111 211 . The reactor fuel is plutonium- 
uranium mixed oxide (MOX) which consists of en- 
riched uranium dioxide UOQ) to 18 w% in 235 U and 
plutonium dioxide PuOQ). Fraction of fissile Pu con- 
tent (( 239 Pu + 241 Pu)/all) is about 16 w% at the inner 
core and about 21 w% at the outer core. 

JOYO reactor operates for 60 days then stops for a 
few weeks in its operational cycle. Therefore, we could 
collect data in both the reactor-on and reactor-off condi- 
tions. The data taken under reactor-off condition were 
used to measure the background. Thermal power of the 
rector was stable at 140 MW during its operation. Fig- 
ure[T]shows time variations of fission rates of main iso- 
topes ( 235 U, 238 U, 239 Pu and 241 Pu) in the fuel. Neutrino 
flux from the reactor core was calculated from available 
measurements ofy6-decay spectra with 2.5 % systematic 
uncertainty U 

3. Neutrino detection principle 

Reactor neutrinos are detected with a liquid scintilla- 
tor formulated from organic oils. Organic oils are abun- 
dant in free protons and the reactor v e react with the 
proton through inverse /?-decay reaction. 



v e + p — > e + n 



(3) 



Figure [2] shows the reactor neutrino flux at JOYO ex- 
perimental site and the cross-section of inverse /3-decay 



152 
153 
154 
155 
156 
157 
158 
159 




Energy ( MeV ) 



Figure 2: Shape of neutrino energy spectrum (arbitrary unit) via in- 
verse y3-decay reactions expected in the detector (solid line). Overlaid 
curves show the reactor neutrino flux (dashed line) and cross-section 
of inverse /3-decay reaction (dotted line). 



reaction together with a shape of the energy spectrum 
in the detector. Number of interactions in the detector 
is determined as a multiplication of the flux, interac- 
tion cross-section and the number of free protons in the 
detector. Energy threshold of the inverse /3-decay in- 
teraction is 1.8 MeV. Cross-section of inverse yS-decay 
reaction is associated with the lifetime of free neutrons 



and calculated precisely with 0.2 % accuracy IU5I1 . 



The detector contains gadolinium-loaded liquid scin- 
tillator (Gd-LS), in which neutrino signals are detected 
by using delayed coincidence technique. A positron ki- 
netic energy and y's from its annihilation are observed 
as the prompt signal. Since the recoil energy of neu- 
tron is small, neutrino energy can be measured from the 
energy of the prompt signal. 



J signal 



1.8MeV + 2w„c 2 



(4) 



Neutrons from inverse /3-decay reactions are captured 
by gadolinium or hydrogen in the Gd-LS mostly after 
thermalization, and y-rays are emitted. Those y-rays 
are detected as delayed signal. In our detector, we ex- 
pect 76.9 % of neutron captures are on 155 Gd or 157 Gd, 
which have more than 10 5 times larger thermal neutron 
capture cross-section than hydrogen II 1611 . y-rays with 
total energy of approximately 8 MeV are emitted from 
a neutron capture on Gd. The mean time difference {At) 
between the prompt and delayed signals is estimated to 
be 46/isec. The background events are strongly sup- 
pressed by requiring coincidence of two signals. 
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Figure 3: Schematic view of experimental setup of the detector. 




164 4. The detector 

165 4.1. Experimental setup 

166 Figure [3] shows a schematic view of the detector. The 

167 detector was constructed at Tohoku University as one 
lee of the R&D programs for the KASKA reactor neutrino 

169 oscillation experiment A 1 Qfl . The detector was moved 

170 to the Joyo reactor site in September 2006 after the 

171 KASKA R&D studies, and was set up on the ground 

172 floor of the reactor building near a delivery entrance, 

173 just at the west outside of the reactor containment ves- 

174 sel. The distance to the Joyo core was 24.3 m. The lo- 

175 cation of the detector is shown in Figure [4] Because 

176 Joyo uses sodium coolant, water was not allowed to be 

177 brought in the building and water shield was not possi- 

178 ble. 

179 The data taking period was from January 2007 until 
iao December 2007. Unfortunatelly the liquid scintillator 
iai deteriorated during the operation and only net 38.9 days 

182 reactor-on data and 18.5 days of reactor-off data were 

183 used for the analysis. The reason of the deterioration is 

184 not clear but we assume the high temperature environ- 

185 ment and N2 bubbling were possible reasons. Details of 

186 detector design and components are shown in the fol- 

187 lowing sections. 

188 4.2. Main detector 

189 The main detector consisted of 0.76 tons of Gd-LS 

190 filled in a transparent acrylic spherical vessel with in- 

191 ner diameter of 1.2 m. The vessel was made of UV 

192 transparent acrylic ACRYLITE(OOO) of MITSUBISHI 

193 RAYON Inc.. Two acrylic hemispheres were made 

194 from 15 mm thick acrylic plates by vacuum form- 

195 ing. The two hemispheres were put together sand- 

196 wiching Viton O-ring at the equator to form a sphere. 



Figure 4: Location of the detector at experimental fast reactor JOYO. 



Parameter 



Value 



199 
200 
201 
202 

203 
204 
205 
206 
207 
208 
209 
210 
211 
212 
213 
214 
215 
216 
217 



Density (20 °C) 
H/C ratio 
Number of Protons (H) 
Light yield 
Gd concentration 
Neutron capture time 



0.838 g/cm 3 

1.94 
6.22x 10 28 
9,400 photon/MeV 
0.05 w% 
46.4yusec 



Table 2: Properties of gadolinium-loaded liquid scintillator used in 
our detector 



There is a 30 cm diameter chimney at the top of the 
sphere. The acrylic sphere was supported by an alu- 
minum stand which stood in a oil pan. The liquid 
scintillator was formulated by diluting the commercial 
Gd- loaded liquid scintillator BC521 (Saint-Gobain) by 
Paraffine oil and Pseudocumene. The compositions of 
the liquid scitillator were, 12.6weight% (w%) Pseu- 
documene (1,2,4-Trimethylbenzene: C9H12), 76.3 w% 
Paraol 850, 11.2w% BC521, and 1.52g/liter of PPO 
(2,5-Diphenyloxazole: C15H11NO) as the fluor. The Gd 
concentration was 0.05 w% (as contained in BC521). 
Paraol 850 is heavy isoparaffin, one of Shell prod- 
ucts. The scintillation light yield was measured to be 
56 % of Anthracene scintillator, which is equivalent to 
9,400 photons/Me V. The Gd-LS was purged by N 2 bub- 
bling with a flowing rate of lOOcc/min during operation 
to reduce the oxygen quenching effect. Properties of 
Gd-LS used in our detector are summarized in Table [2] 
The scintillation lights from Gd-LS were measured 
by 16Hamamatsu R5912 8-inch photomultiplier tubes 
(PMTs) mounted on the surface of the acrylic vessel. 




Figure 5: Picture of sphere shape acrylic vessel equipped with 16 
PMTs on the surface. Gd-LS is not filled at the time of this picture. 

2ie Each PMT was covered by a mu-metal skirt which was 

219 used for Kamiokande PMT long time ago. Figure [5] 

220 shows a picture of the acrylic vessel with 16 PMTs 

221 on the surface. The PMT was put in a acrylic hous- 

222 ing cylinder and the space between PMT surface and 

223 acrylic sphere were filled with RTV rubbers (Shin-Etsu 

224 Silicones KE103, KE1052). The photo-cathode cover- 

225 age was approximately 10 % . 

226 4.3. Cosmic-ray veto counter and detector shielding 

227 Since the detector was set up at ground level, cosmic - 

228 ray muon flux was large. In order to reduce the 

229 cosmic -ray muon background, the main detector was 

230 surrounded by a veto counter system. This system 

231 consisted of two layers of 1 cm thick plastic scintilla- 

232 tor plates equipped with wavelength shifter and PMT 

233 for the readout. Top of the detector and the north and 

234 south sides were fully covered by scintillator layers, 

235 while only the upper half were covered for the east and 

236 west sides. Veto efficiency of cosmic -ray muons by the 

237 counter system was estimated to be 92 % from MC sim- 
23e ulation including the acceptance. The cosmic -ray veto 

239 signal rate was about 2 kHz. 5 cm thick lead blocks cov- 

240 ered the bottom area of the oil pan and 6 mm thick lead 

241 sheets backed by wooden boards covered the side of the 

242 detector housing. 5 cm thick paraffin blocks were ar- 

243 ranged outside of the detector to suppress fast neutrons 

244 induced by cosmic muons. Figure [6] shows a picture of 




Figure 6: A picture taken at the experimental site. The detector is in 
the black cottage at the left hand side. The reactor containment vessel 
is behind the concrete wall at the right hand side. 

245 the detector at the site. The size of the detector was 

246 roughly 2.5 mx2.5 mx2 m(H). 

247 4.4. Data acquisition system 

248 CAMAC standard electronics modules were used for 

249 data taking. Figure [7] shows the schematic view of the 

250 data acquisition system. Signal from each PMT was 

251 divided into four. The first one was fed into ADC, 

252 by which integrated charge within 200nsec time win- 

253 dow was measured. The second signal was used to 

254 make a common trigger for ADC and TDC, which was 

255 made from a discriminated analog sum of all PMT sig- 

256 nals. The threshold for the common trigger was set at 

257 3.5 MeV. If cosmic -ray veto counter had a signal within 

258 100 /vsec before the common trigger was made, the trig- 

259 ger was canceled. The third signal provided a stop sig- 

260 nal to TDC, which measured the timing of the PMT 

261 hits. The last signal was fed into another ADC for pulse 

262 shape discrimination (PSD) study aiming to identify fast 

263 neutron background, although the PSD was not used for 

264 the study described in this paper. In order to collect 

265 delayed coincidence signals from neutrino interactions, 

266 lower trigger threshold at 2.5 MeV was applied to the 

267 delayed signals for 100 //sec after a prompt trigger was 

268 created. Time interval between the first and second trig- 

269 gers was measured by counting a 100 MHz clock signal 

270 by a CAMAC scaler and the data were saved along with 

271 the ADC and TDC data for each trigger. If a trigger for 

272 delayed signal was not generated within 100 //sec, the 

273 data acquisition system was back to the normal mode 

274 with 3.5 MeV threshold. In addition to this delayed trig- 

275 ger, we also took data with single trigger at 0.6 MeV 

276 threshold for the background study. During the data tak- 

277 ing at JOYO fast reactor site, the single trigger rate was 

278 about 300 Hz. Mean dead time of the data taking inclu- 

279 sive of the cosmic muon veto time was 38 %. 
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The readout and monitoring system in this experi- 
ment needed to be simplified due to limited access to 
the experimental area. Therefore, we constructed our 
DAQ software system in a CAMAC CC/NET 01 to 
read the data from ADC and TDC modules and used a 
trigger system installaed into a NIM FPGA module. 

As JOYO is a fast reactor and uses sodium as mod- 
erator, the experimental area also needed to be kept off 
water and high humidity. In addition, the detector us- 
ing liquid scintillator, which generates organic gas, was 
placed in a large box sealed with black vinyl sheets. In 
order to keep safety of the experimental area during the 
operation of the experiment, we constructed a monitor- 
ing system in a Linux computer and kept watching tem- 
perature, humidity and density of oxygen and organic 
gases. The experiment and monitoring data were auto- 
matically sent to a 220 km distance remote site, Tohoku 
Gakuin University. We built a secure network on the 
internet by IPsec VPN architecture over IPv4 protocol, 
which enabled an experiment shift person to check the 
condition of DAQ and experiment area remotely [18]. 



4.5. Monte Carlo simulation 

The collected data were compared with the Monte 
Carlo (MC) simulation based on Geant4 (version 
4.9.0.pl). Geant4 is a toolkit which provides a 
calculation of particle tracking in materials fl^l . 
For the hadronic interaction process, QGSP_BIC_HP 
model i20Tl was employed in Geant4. It comprehends 
from low energy region under 20 Me V such as behavior 
of thermal and fast neutron to high energy region such 
as interactions between cosmic-ray muons and materi- 
als around it. Trajectory of optical photons emitted in 
the Gd-LS was simulated considering the optical pro- 
cess including attenuation and scattering. 

Corrections for PMT responses and energy calibra- 
tion were carried out by putting a 60 Co y-ray source 
inside the detector. Energy was reconstructed from 
the total observed charge by 16 PMTs in which cor- 
rection to the acceptance and attenuation length in the 
liquid scintillator were taken into account. Figure [8] 
shows a reconstructed energy spectrum from the data 
taken with a 60 Co y-ray source at the detector center. 
The 60 Co source mainly emits two gamma rays with 
1.17MeV and 1.33 MeV energies. A large peak in Fig- 
ure[8]is made from the gamma rays with 2.5 MeV total 
energy. The energy resolution estimated from the peak 
at 2.5 MeV was 20 %/ V£(MeV). 

In addition, the measured data with 24I Am- 9 Be (a, 
n) neutron source at the detector center were used to 
tune the quenching effects of protons recoiled by neu- 
trons parametrized by Birks' constant kg I2H1 and eval- 
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Figure 7: The schematic view of the data acquisition system. 
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Figure 8: Reconstructed energy spectrum from data taken with Co 
y-ray source at the detector center. 
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Figure 9: Comparison between reconstructed energy spectra of the 
prompt trigger events for the reactor-on and off. Black and gray his- 
tograms show the observed data for a day live-time under reactor-on 
and off conditions, respectively. There is an excess around 8 MeV of 
the distributions attributed to thermal neutron capture on Gd. 

331 uate the neutrino MC simulation. The Birks' constant 

332 of our Gd-LS was estimated to be 0.07 mm/MeV from a 

333 comparison of the measured energy spectrum to the MC 

334 simulation. 

335 Not only the neutrino signal events, but also various 

336 background events were generated by the MC simula- 

337 tion and compared with the observed data. Those back- 

338 ground events included cosmic-ray muons and the muon 

339 decay, fast neutrons and environmental y-rays from de- 

340 cay chains of 238 U and 232 Th series and 4() K decays. In 

341 addition to the fast neutron and environmental y-rays 

342 generated inside the detector, those from outside of the 

343 detector were also considered in the MC simulation. 

344 5. Measurement of background spectrum 

345 Major background sources in this experiment were 

346 environmental y-rays and cosmic -ray muons. The en- 

347 vironmental y-rays are emitted by radioactive isotopes 

348 contaminated in the detector and materials around the 

349 detector. These y-rays are produced through the de- 

350 cay chains of 238 U and 232 Th series, and decay of 4() K. 

351 The energy of y-rays ranges up to 2.6 MeV. However, 

352 there were y-ray contaminations above the discrimina- 

353 tor threshold level of 3.5 MeV due to the energy res- 

354 olution tail. The main source of the y-rays was con- 

355 sidered to be concrete walls surrounding the detector. 

356 Cosmic-ray muons have wide energy range over GeV 

357 scale. High energy muons generate fast neutrons by in- 

358 teractions in materials composing the experimental site 

359 and fast neutrons turned out to be the severest back- 

360 ground for neutrino signals in this experiment. Most 



361 of the background events produced by muons were ex- 

362 eluded by the delayed coincidence technique but there 

363 were still remaining backgrounds even after requiring 

364 it. Those background events were further reduced by 

365 the data analysis as explained in later sections. 

366 Huge number of neutrons were produced in the re- 

367 actor core, and a very small fraction of them could 

368 reach the neutrino detector passing through materials 

369 constructing the Joyo or crevices in the materials. Low 

370 energy neutrons are detected as gamma-rays emitted via 

371 neutron capture on Gd. Figure|9]shows a comparison of 

372 the reconstructed energy spectra for the reactor-on and 

373 off. An excess in reactor-on was found around 8 MeV, 

374 which was considered to be made by the low energy 

375 neutrons from the core. The excess rate between 4 MeV 

376 to 14 MeV was 6.1 Hz. Because energy distribution of 

377 the thermal neutron was unknown, it was impossible 

378 to precisely calculate the detection efficiency. There- 

379 fore we carried out rough estimation of the thermal neu- 

380 trino flux. Assuming the detection efficiency of 8 MeV 

381 gamma rays was 50% and considering naively the total 

382 cross section for the neutrons to the detector was equiva- 

383 lent to the surface area of the acrylic sphere, the neutron 

384 flux at the detector site can be estimated to be approxi- 

385 mately 10~ 4 /cm 2 /sec. This neutron flux is too small to 

386 be detected by usual neutron counters. 

387 The background events which satisfy the delayed co- 

388 incidence condition were classified into two categories, 

389 namely accidental and correlated backgrounds. The ac- 

390 cidental background consists of two independent back- 

391 ground events which accidentally occur within the de- 

392 layed coincidence time window. The main source of 

393 such background events are environmental y-rays fol- 

394 lowed by cosmic -ray muons. The correlated back- 

395 grounds are caused by a continuous physics process. 

396 Those processes include decays of cosmic -ray muons 

397 inside the detector and fast neutrons induced by cosmic- 

398 ray muon followed by neutron capture on Gd in the 

399 detector. In the former case, cosmic-ray muon causes 

400 background to the prompt signal, and an electron from 

401 the muon decay (Michel electron) is identified as the de- 

402 layed signal. In the latter case, recoil protons caused by 

403 a fast neutron are detected as a prompt signal, and the 

404 y-rays from neutron capture on Gd are identified as a 

405 delayed signal. Some of radioactive isotopes produced 

406 by cosmic -ray muon interactions in the detector cause 

407 coincidence signals in its decay chain and can be con- 
4oa sidered as background source |22|. However, the pro- 

409 duction rate of such isotopes was negligibly small com- 

410 pared to the other background sources in this experi- 

411 ment. 

412 Figure \W\ shows the reconstructed energy spectrum 
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Figure 10: Reconstructed energy spectra of the prompt trigger events above ().6MeV of the threshold level with different energy ranges. Points 
show the observed data for 309 sec live-time taken in reactor-off condition. Overlaid histograms show the expected neutrino signal and background 
energy spectrum with the contributions from each background source. 



413 of the observed data with 0.6MeV of the threshold level 

414 together with the expected reactor neutrino signal and 

415 background events from the MC simulation, in which 

416 delayed coincidence cut condition was not required yet. 

417 Simulation of cosmic-ray muon background was based 

418 on flux measured in [23], and correction factor 0.72 was 

419 applied from a fit to the data for the high energy re- 

420 gion between 20 MeV and 140 MeV. Fast neutron flux 

421 was obtained as 0.63 neutrons/cm 2 /sec above 10 MeV 

422 (equivalent to visible electron energy of 3 MeV in the 

423 detector) from a fit to the distribution of time interval 

424 (Af) between the prompt and delayed signals. As is 

425 shown in Figure [10l environmental y-ray background 

426 is dominant for the energy below 6 MeV before the de- 

427 layed coincidence condition is applied. We assumed 

428 typical concentration of radioactive isotopes in concrete 

429 materials around the detector, 2. 1 ppm of 238 U, 5.1 ppm 

430 of 232 Th and 1.4 ppm of 40 K [24]. The observed en- 

431 ergy spectrum was in reasonable agreement with the 

432 expected background spectrum from the MC simula- 

433 tion after normalization corrections were applied to the 

434 cosmic -ray muon and fast neutron flux. As is shown 

435 in this figure, the background level is 10 s ~ 10 7 times 

436 higher than the neutrino signals before the delayed co- 

437 incidence. In order to measure reactor neutrino sig- 

438 nals significantly over the background, reduction of the 

439 background events by the delayed coincidence tech- 

440 nique and further selections are necessary. 

441 6. Neutrino event selection and results 

442 We used the data taken with delayed coincidence trig- 

443 ger to search for neutrino events. In order to reduce the 

444 backgrounds, following selection criteria were applied 

445 to the data and the MC simulation. 



446 • 4.5 < E prompt - 7 MeV and 

4.5 < Edelayed < 1 1 MeV, 

448 where E prompt and E^eiayed are the reconstructed energy 

449 for the prompt and delayed signals, respectively. Fig- 

450 ure QT] shows the reconstructed energy spectra for the 

451 prompt and delayed signals. The lower cut value at 

452 4.5 MeV was set to reject environmental y-rays, while 

453 the higher cut values were set at 7 MeV and 1 1 MeV to 

454 select prompt signal shown in Figure [2] and total 8 MeV 

455 y-rays from neutron capture on Gd, respectively. Es- 

456 pecially, the delayed energy selection is very effective 

457 for rejection of Michel electron events with 53MeV of 

458 end point, in contrast to the fast neutron events induced 

459 by cosmic muons which have a similar distribution to 

460 neutrino delayed signals. 

461 • 2.5 < Af < 60/Ksec 

462 Figure[12]shows the time difference between the prompt 

463 and delayed signals (Af) for the events after the energy 

464 cuts. Michel electron events have distribution follow- 

465 ing 2.2//S of muon life time, while the fast neutron and 

466 neutrino events have corresponding distributions with 

467 the decay time of 46.4yus, which is determined by Gd 

468 concentration in liquid scintillator. Although the coin- 

469 cidence condition within 100/ks time window was re- 

470 quired in the data acquisition, further cut was applied in 

471 the analysis. The lower limit was used to reject remain- 

472 ing Michel electron events after the energy cuts, while 

473 the upper limit was set to collect enough neutrino events 

474 with 46.4/is decay time. 

In order to further reduce background events re- 
mained after the energy and Af cuts were applied, 
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Figure 11: Reconstructed energy spectra for the prompt (left) and delayed (right) signal candidate events. Points show the observed data for 7.4 
hours live-time under reactor-off condition. Overlaid histograms show the expected spectra of neutrino signals and background events from MC 
simulation. 



charge balance (CB) is defined as follows: 



483 
484 
485 



CB 



16 



(si 6 , (er) 2 ) 



(5) 



where Q c j orr is observed charge from z'-th PMT after gain 
correction was applied. This variable becomes large for 
the external background events, such as environmental 
gamma rays and fast neutrons, with the vertex position 
close to the surface of the detector while it is smaller 
for events occuring near the center. So CB cut corre- 
sponds to a kind of fiducial volume cut. However when 
the vertex position is too close to the surface and be- 
tween the near PMT surfaces, the CB becomes smaller 
and the vertex position mimics a place around the cen- 
ter of the detector, because solid angles to the near PMT 
surfaces from the vertex position are narrow and number 
of photoelectrons for the near PMTs becomes less. Fig- 
ure [13] shows the CB distributions of reactor neutrino 
and background events after energy and At cuts were 
applied. Distributions for the background events have a 
valley at CB ~ 1.3. To maximize signal over noise ratio 
(S/N), cut conditions for CB were defined as: 



0.8 < CB 



prompt 



< 1.4 and 0.8 < CB delaxed < 1.4 



Even after the energy, Af and charge balance cuts 
were applied, the remaining background events were 
still hundred times larger than the reactor neutrino 
events. The dominant component of the remaining 
background events was fast neutrons generated from 
muons in materials around the detector. Figure [14] 




' 10 20 30 40 50 60 70 80 90 100 

At (us) 



Figure 12: At distributions for events which satisfy the energy cut: 
4.5 < Eprampi < 7 MeV and 4.5 < E delayed < 11 MeV. Points show 
the observed data for 7.4 hours live-time taken in reactor-off condi- 
tion. Overlaid histograms show the expectations of neutrino signal 
and background events from MC simulation. 



Selection Criterion Event rate (/day) S/N ratio 



Trigger level 


162 




Energy cut 


8.89 


1/1197 


Coincidence cut 


7.01 


1/1009 


Charge balance cut 


0.988 


1/128 


Vertex <p cut 


0.494 


1/34.6 



Table 3: Effect of selection criteria on the reactor neutrino events in 
the detector. Each cut condition is described in Section[6] 
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Figure 13: Charge balance (CB) distributions for the events which satisfy energy and At cuts. Left-hand and right-hand figures show the distributions 
for prompt and delayed signals, respectively. Points show the observed data for 7.4 hours live-time taken in reactor-off condition. Overlaid 
histograms show expected distributions of neutrino signal and background events from MC simulation. 



502 
503 
504 
505 
506 
507 
508 
509 
510 
511 
512 
513 
514 

515 
516 



shows the vertex (p distributions for the prompt and de- 
layed signals after energy, At and charge balance cuts 
were applied. <p is the azimuthal angle in spherical polar 
coordinates as the z-axis vertical to the ground. Then 
<p — was defined as north of the detector. The vertex 
position is reconstructed by a fit with expected charge 
of each PMT based on the scintillation light yield and 
a solid angle to the PMT from the vertex position. It 
is expected that the vertex positions of the neutrino in- 
teractions distribute uniformly in the detector. On the 
other hand, the vertex tp distributions of the data are not 
flat due to asymmetric arrangement of building materi- 
als and paraffin shields. Therefore, we applied the fol- 
lowing cut based on the vertex <p position to maximize 
S/N ratio: 

• -100° < <p prompt < 100° and 
-100° < 4> delayed < 100° 

The neutrino event selection cuts and expected event 
rates are summarized in Table [3] 

6.1. Result of the reactor neutrino event selection 

The result of the neutrino event selection is summa- 
rized in Table |4] The accidental background event rates 
were estimated from the measurements of single back- 
ground event rates and subtracted from the total event 
rates as shown in Table [4] Then, the difference be- 
tween the event rates for reactor-on and reactor-off was 
calculated. Errors in Table |4] are only the statistical 
ones. After reactor neutrino event selection and subtrac- 
tion of accidental background and reactor-off data, event 
rate for the neutrino candidate events from 38.9 days of 



Reactor-on Reactor-off A(on - off) 



Live-time 

Total 
Accidental 
Correlated 



38.9 days 
19.0±0.7 
2.34±0.24 
16.7±0.7 



18.5 days 
17.1 + 1.0 
1.52+0.29 
15.6+1.0 



1.93+1.19 
0.82+0.38 
1.1 1 + 1.24 



Reactor v MC 



0.494 



Table 4: Observed event rates (events/day) and the statistical uncer- 
tainties after the neutrino event selection criteria were applied. Acci- 
dental background event rates were estimated by single background 
event rate. Correlated event rates were obtained by subtracting the 
accidental BG event rate from the total event rate. 



reactor-on data and 18.5 days of rector-off data was ob- 
tained as 1.1 l + 1.24(stat.)+0.46(syst.) events/day, while 
the expected neutrino signal event rate from the MC 
simulation was 0.49 events/day. The systematic uncer- 
tainties were estimated considering the uncertainties in 
energy resolution, energy scale and vertex reconstruc- 
tion. Figure Q3] shows the prompt energy spectrum af- 
ter all selections except for the prompt energy cut. The 
measured excess rate was consistent with the expected 
neutrino rate from the MC simulation, but also consis- 
tent with zero within the systematic error. So the obser- 
vation of neutrinos from experimental fast reactor JOYO 
has not been statistically established in this measure- 
ment. A new design of the detector is described in the 
following section, in which sensitivity to the measure- 
ment of reactor neutrinos was estimated based on the 
observed background event rates shown in this paper. 
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Figure 14: Reconstructed vertex <p distribution after the cuts of energy, At and CB were applied. Points shows the data for 18.5 days live-time taken 
in reactor-off condition. An overlaid histogram shows the expectation of neutrino signals. 
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Figure 15: Reconstructed energy spectrum after neutrino event selec- 
tion criteria except for the prompt energy cut was applied. Points show 
the data with the statistical errors in which energy spectrum measured 
for reactor-off is subtracted from that taken for reactor-on. Accidental 
background event rates were estimeted from the data and subtracted. 
Boxes show the expected reactor neutrino energy spectrum from MC 
simulation with the MC statistical errors. 

547 7. New detector design for the next experiment 

548 Problems found in the JOYO experiment were fol- 

549 lowing. 

550 (1) Statistics of the data was limited because the long 

551 term data taking was impossible due to degradation 

552 ofGd-LS. 

553 (2) Energy threshold level could not be sufficiently 

554 lowered due to large amount of environmental y- 

555 rays from outside of the detector. Therefore, the 
sse neutrino detection efficiency was obliged to be low. 



557 (3) Performance of event vertex reconstruction was 

558 not enough to distinguish external background 

559 events entering the detector due to the same prob- 

560 lem caused by the detector structure as in the de- 

561 scription for CB in Section|6] 

562 (4) Fast neutron background level was too high against 

563 the neutrino signals, for which S/N ratio was 

564 0.029 . 

ses We designed a new detector for the next experiment 

see solving these problems, and estimated the sensitivity to 

567 reactor neutrino measurements based on the MC simu- 

sea lation. 

569 We are considering a liquid scintillator with high 

570 Pseudocumene concentration above 99 w% as possible 

571 candidate to solve the Gd-LS degradation problem. The 

572 high aromatic concentration is supposed to stabilize the 

573 Gd-LS. Experimental studies of the long-term stabil- 

574 ity and characteristics of the proposed Pseudocumene 

575 based Gd-LS are necessary. In addition, the detector de- 

576 sign needs to be improved to suppress fast neutron back- 

577 ground. The new detector will consist of two concen- 

578 trie sphere vessels. The inner vessel contains the Gd-LS 

579 as target of reactor neutrinos. The outer vessel is filled 

580 with paraffin oil with no scintillation light emission, and 

581 works as shield against fast neutrons. In the MC simu- 

582 lation, we assumed the same target volume with Gd-LS, 

583 surrounded by 20 cm layer of buffer oil. Environmen- 

584 tal y-rays are also reduced by the outer layer, by which 
sss we estimate that the energy threshold can be lowered to 
sse 3 MeV with the same trigger rate as the measurement at 
587 JOYO. The scintillation lights are viewed by 24 10-inch 
ssa PMTs isotropically arranged on the surface of the outer 
589 vessel, providing 1 1 % photo-cathode coverage close to 
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590 JOYO detector. The PMT surfaces are away from the 

591 target vessel with the interval buffer region by which 

592 performance of vertex position reconstruction can be 

593 improved especially for those close to the surface of the 

594 target vessel. The reconstructed vertex radii in the polar 

595 coordinate system are used for rejection of the external 

596 background events. 

597 According to the study using the MC simulation, we 

598 expect the S/N ratio is improved from 0.029 to 0.093, 

599 with about 16 times larger neutrino selection efficiency 

600 by the new detector design. If we put the same detector 

601 at MONJU rector site M25I1 . of which thermal power is 

602 approximately five times higher than JOYO reactor, the 

603 S/N ratio is further improved to 0.48 with 41 events/day 

604 of neutrino observations. Significance of fast reactor 

605 neutrino observation reaches 2 standard deviation after 

606 12 days of reactor-on and off live-times at JOYO reactor 

607 site or 1 day at MONJU reactor site. 



608 8. Conclusions 



609 We carried out an experimental study of fast reac- 

610 tor neutrino detection and measured the background 
en spectrum at fast reactor JOYO using a compact detec- 

612 tor. The observed reactor neutrino candidate signal 

613 was l.ll+1.24(stat.)±0.46(syst.)events/day after sub- 

614 traction of background events while the expected re- 
els actor neutrino event rate from the MC simulation was 

616 0.49 events/day. As a result, the first observation of fast 

617 reactor neutrinos was not statistically established in this 

618 measurement. On the other hand, various background 

619 sources at the ground level nearby the reactor were stud- 

620 ied in detail and those backgrounds were found to be 

621 reproduced well by our MC simulation. These back- 

622 ground studies will be useful not only for the R&D of 

623 future reactor neutrino oscillation experiments but also 

624 for the development of compact reactor neutrino detec- 

625 tor as a remote monitor for plutonium breeding. 

626 A design concept of a new detector and its sensitivity 

627 to the observation of fast reactor neutrinos were also de- 

628 scribed in this paper based on the knowledge acquired 

629 by the measurement at JOYO. We expect the S/N ratio 

630 will be improved from 0.029 to 0.093 by the new detec- 

631 tor design, and it is further improved to 0.48 if we put 

632 the detector at MONJU reactor site at the same distance 

633 from the core. Expected reactor neutrino signal by the 

634 new detector design is 8.0 events/day and 41 events/day 

635 at JOYO and MONJU reactor site, respectively. 
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